Optical spatial solitons in soft-matter: mode coupling theory approach 
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We predict that spatial self-trapping of light can occur in soft matter encompassing a wide class 
of new materials such as colloids, foams, gels, fractal aggregates etc. We develop a general nonlocal 
theory that allows to relate the properties of the trapped state of Maxwell equations to the measur- 
able static structure factor of the specific material. We give numerical evidence for stable trapping 
in fractal aggregates and suggest also the possibility of soliton spectroscopy of soft-matter. 



Self-trapping of light beams, predicted forty years ago 
Q, is still a subject of great interest 0, Q Obser- 
vation of optical spatial solitons (OSS) at low (down 
to mW) power levels, demonstrated in photorefractive 
or liquid crystals, makes OSS attractive candidates in 
several applications of emerging photonics technology, 
0,0,0; and have driven successful efforts to understand- 
ing the role played by specific material properties such 
as anisotropy and/or nonlocality in self-trapping mecha- 
nisms. Yet, the description of light trapping seems still 
strongly material dependent, and, as such, cannot be ap- 
plied to a whole class of condensed matter, namely soft 
matter (SM), which encompasses the important case of 
bio- matter where OSS can find new applications (e.g., 
laser surgery, optical manipulation of nano-particlcs). 
Roughly speaking, softness is generally due to a charac- 
teristic mcsoscopic (i.e., larger than atomic) length scale 
of the constituents, and an energy scale comparable to 
room temperature thermal energy. As a consequence, 
SM properties can be easily tailored via external field of 
different (mechanical, electrical, magnetic, thermal, ...) 
origin. SM includes colloidal suspensions, emulsions, and 
foams (all involving different constituents in a host fluid) , 
and typical examples are polymers in a liquid, glues, li- 
posomes, blood, and all sort of bio-matter [jj. 

Historically, the use of aerosols and water suspension 
of dielectric spheres as nonlinear media dates back to 
the early 80's 0, Q. However, in these materials elec- 
trostriction has been described in the framework of sim- 
ple models that treat the diluted constituents as a gas of 
non-interacting particles, in terms of local index change 
of the Kerr type (An = ri2l where / is the local intensity) 
IH0E1 As a consequence existence of stable OSS in two 
transverse dimensions is ruled out by the occurrence of 
the well-known catastrophic self- focusing instability ^3 • 
Phcnomenologically, stabilization can be expected from 
the index saturation arising from the maximum pack- 
ing fraction of the dispersed particles. However, before 
reaching such a regime, the physics of these materials, 
and in general of other SM systems, is well-known to 
be affected by particle-particle interactions. In particu- 
lar, this occurs whenever the particle-particle correlation 
function g(r) is structured on a length scale comparable 



to the laser beam waist. Under these conditions a non- 
local model for self-focusing in SM must be considered. 

In this letter, we propose a novel general approach to 
stationary self-focusing, linking for the first time the elec- 
trostrictive nonlocal response of SM to its static structure 
factor S(q) (roughly speaking, the Fourier transform of 
the particle-particle correlation function), usually mea- 
sured by means of scattering experiments. This allows 
us (i) to predict stable propagation of two-dimensional 
OSS in a new wide class of condensed matter, and (ii) to 
assess the importance that ultra-focused laser light can 
have to investigate the properties of SM. Both issues are 
of paramount importance in order to go towards a more 
general description of solitons in complex media (fractal 
aggregates, structured and supercooled liquid, etc.) and 
their application in bio-photonics as well as to develop a 
new spectroscopic tool for the investigation of SM prop- 
erties. 

Assuming a linearly polarized beam and exploiting 
isotropy of system, we start from the unidirectional scalar 
wave equation written for a monochromatic beam with 
complex amplitude £(x,y,z) propagating along z \\\\ 
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where k 2 = uj 2 e(bj)/c 2 = ijp'n^jc 2 , no is the SM bulk re- 
fractive index, and = d 2 + d 2 . We further assume 
that the nonlinear polarization is responsible for a refrac- 
tive index change An, i. e. P n i{£) = A^f = 2noAn£ , 
which is dominated by an electrostrictive contribution 
(fast electronic nonlinearities as well as index change due 
to thermal heating by optical absorption are usually neg- 
ligible, and we also neglect scattering). This legitimates 
our scalar (polarization-indipendent) approach [l2T |. re- 
gardless of beam spectral content [at variance with Kerr 
effect of electronic origin where the non-paraxial regime 
at very high intensities requires to account for vectorial 
effects (see e.g. fl^)]. 

Since the electrostriction An = p(dn/dp) Po is propor- 
tional to the particle number density change p (from 
equilibrium value po), Eq. (JJJ must be coupled to an 
evolution equation for p. A widely accepted and largely 
applicable theory for SM is the so-called Mode-Coupling 
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Theory (MCT) 003113 which relics on the so-called 
Zwangig-Mori formalism (given some observable, like p, 
it allows to write closed equations for it and for its cor- 
relations functions By exploiting MCT (see also 
|l8l|). we find that the density perturbation p = p(q,t) 
in Fourier space q = {q x ,<ly,<lz), obeys the dynamical 
equation: 



/(?,*) 
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q 2 £m(t-t')'p(q,t')dt' 
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where tilde denotes 3D spatial Fourier transform, rj = 
(/Lto/eo rl o) 1 ^ 2 i s the impedance, ks is the Boltzman con- 
stant, T is the temperature, m(t) is the memory function 
of the system (for a simple liquid m(t) is a Dirac delta 
function times the viscosity), and S(q) is the static struc- 
ture factor. In Eq. @ / is a Langevin term describing 
random forces (see e.g. (l8j1- 

Eq. J5J) has been always considered without the de- 
terministic forcing term wei ghte d by the electrostrictive 
coefficient 7 e = p(de/dp) Po , and provides one the 
most successful approaches to structural phase transi- 
tions of soft-matter. Here we extend it by accounting for 
the presence of an external optical field [coupled through 
Eq. JJJ] , which induces an electrostrictive force with po- 
tential proportional (in configurational space) to V 2 /. 
This model can be thou ght of as a generalization of the 
acoustic wave equation [12j , which has been previously 
employed to determine the electrostrictive correction to 
electronic nonlinearity of silica glass ^| • The latter case 
is retrieved in our approach for S(q) = constant and a 
memory-less response [i.e. m{t) oc 8{t)\. In the general 
case considered here, MCT accounts for the elastic defor- 
mation of a medium made of interacting particles which 
results in a non-homogeneous response weighted by S(q). 
Additionally, though here we deal with electrostriction, 
MCT is a powerful approach that can be generalized to 
account for other types of nonlinearity, e.g. reorenta- 
tional mechanisms, by looking at different observables. 
As such, it provides a general framework for studying 
nonlinear optics in SM, going beyond the idealized Kerr 
(local) limit H Q (which is, nevertheless, correctly re- 
trieved for non- interacting particles, as shown below). In 
the following, we address specifically the properties of 
spatial solitons. 

Assuming that the random fluctuations arc negligible 
with respect to the driving electrostrictive term, the sta- 
tionary state solution (dt = 0) of Eq. J2J) yields 



(3) 



which shows that S(q) plays the role of a transfer 
function from the optical intensity to the density. 
Incidentally, Eq. J2J can be also obtained by starting 
from a different model of SM employing so-called 
generalized hydrodynamics equations |2(| such as those 



typically adopted to modelling inelastic light scattering 
spectra (ISTS) |21[. Obviously, in real space, Eq. Q 
corresponds to a differential equation for p = p(x,y,z). 
While this equation, as well as Eq. QJ, involves three 
dimensions, OSS imply by definition a z-independent 
intensity. Therefore, in this case, p does not depend on 
z and Eq. coupled to Eq. QJ, wm be interpreted 
henceforth as a 2D transverse equation with q = (q x , q v ) 



Equations (lll.'ill allow us to develop a general non- 
local model for trapping in SM. First we consider the 
paraxial (Fresnel) regime, which corresponds to expand- 
ing the transverse operator in the g-Fourier transform 
JU as 



of Eq. 



fc[l - q 2 /(2k 2 )]. Equation 



t|T~|) becomes, in terms of the slowly varying envelope 
E(x,y,z) = £(x,y,z)exp(-ikz), 
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In Eq. Q the nonlinearity arises from the term pE. In- 
deed, back-transforming Eq. @ to real space for p and 
assuming azimuthal symmetry, we obtain, after some al- 
gebra, a self-consistent nonlinear non-local wave equation 



dE f°° 
i2k— + V 2 ± E + X E G(r,r')\E(r',z)\ 2 r'dr' 
Jo 



0, (5) 



where r = \J x 2 + y 2 , and we have defined the kernel 

G(r,r') = ^ ^lj (Q r )MQr') QdQ, (6) 
Jo 

and the coefficient x = k 2 (dn/dp) Po 'y e So/2kBTno (note 
that S(q) is scaled to So = 5(0) representing the ratio 
between the compressibility of the material and that of 
the ideal gas [12] )■ We seek for bound states of Eqs. QEj) 
in the form E(z,r) = (xwo)~ 1 ^ 2 u(a)exp(i(3Q, where /? 
is the nonlinear correction to the wavevector k, which is 
determined sclf-consistenly in the numerical simulations, 
and a = r/wo, C = z / z o = zj2kw\ are dimensionless 
radial and longitudinal variables, in units of beam width 
Wo and diffraction (Rayleigh) length zq, respectively. The 
OSS (bound state) profile u(a) obeys the non-local eigen- 
value equation (we set Wq^ 7 ! = = d 2 / da 2 +a~ 1 d/ da) 



0u + u / g(cr, a')u 2 (a / ) a 1 da' = 0, 
Jo 



(7) 



where the kernel g(a,a') can be obtained (at least nu- 
merically), once 5(Q) is known, from the integral 



9(cr,a')= r 
Jo 



S(6/w ) 
S 



J (a'6)J {a9) Odd. (8) 



From Eqs. (|7I8|I the ideal local (Kerr) limit is recovered 
for 5(g) = So = constant, which yields g(a,a') = 5{a — 
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cr')/cr, and hence a pure Kerr law with a nonlinear index 
that, in the case of the dielectric spheres, turns out to be 



n 2 



PqSq 
2e h J k B T' 



(9) 



This limit, however, is well known to lead to unstable 
(so-called Townes after Ref. Q) OSS. Conversely, in the 
general case S(q) ^ constant, we expect solutions of Eqs . 
((7151) to be stabilized by non-locality [2II2I El ElE^ . 
We are also naturally brought to consider deviation from 
paraxiality due to strong focusing, and argue for the ex- 
istence of OSS in this case. In fact, narrow OSS in SM 
may be important both in specific applications (e.g. laser 
surgery) and in order to establish OSS as a mean for prob- 
ing the static structure factor S(q) of SM when the lat- 
ter extends to high spatial frequencies. Deviations from 
paraxiality can be accounted for by considering the next 
order in the expansion of the transverse operator in Eq. 
(|T|). By adopting the normalization employed for the 
paraxial case, we cast the new bound state equation in 
the form 



V>- £ V>-/?u- 



i / g(<T,<r')u 2 (a') cr'da' = 0, (10) 
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where the degree of non-paraxialiy is measured by a sin- 
gle dimcnsionlcss parameter e = (A/ 'Airnw^) 2 fixed by the 
ratio between wavelength A and beam width scale Wq. 

In order to discuss various OSS supported by different 
types of SM, we make specific examples. We start consid- 
ering hard spheres in a host liquid (solvent). In the limit 
of diluted, non interacting spheres, S{q) is constant and 
this yields, once again, unstable OSS. In a more refined 
approximation, S(q) can be described by a parabolic law 
in the framework of the Percus-Ycvick model 221 



S(q) = S + Kq 2 



After Eq. (J3J, the corresponding expression for p (as 
stated, we assume p to follow adiabatically / along z) is 



p(r, z) 



2k B T 



[S I(r,z)-KVll(r,zj\ 



[12) 



which, once inserted in Eq. Q), gives a model for 
weakly non-local solitons that has interdisciplinary inter- 
est (plasma physics, matter waves, transport in DNA, see 
Ref. [2^ and references therein) . Stable soliton solutions 
of this model have been reported [23| and, in this context, 
represent 1+2D OSS in SM, when its static structure fac- 
tor can be well approximated by Eq. (|TT)l . To this end, 
consider that, using Percus-Yevick model, the parabolic 
approximation of S(q) breaks down around qr s = 5. 
Since q can be reasonably estimated to be q ~ w^ 1 , the 
weakly non-local model starts to loose its validity when 
the spheres have size comparable with beam width wo (in 
this regime, a microscopic description of the molecular 



dynamics is needed). Viceversa, when Eqs. i|llll2ll hold 
valid, the nonlocality that, generally speaking, provides 
the stabilizing mechanism of OSS against catastrophic 
sclf-focusing |23j stems from the particle-particle corre- 
lation function <?(r), which is proportional to the Fourier 
transform of S(q) — 1 j^, as anticipated. Importantly, 
since S(q) is not uniform, stable OSS not only exist, but 
provide information on the material [for a given optical 
power, the width of OSS is determined by the constant 
K in Eq. fTTfl. 

Further models for S(q) can be discussed. A very in- 
triguing case is that in which the suspended particles of 
colloidal SM develop self-similar aggregates with fractal 
dimension D, described by the function [2^ 



S(ff) = l + 2 ^( 



sin[(D - 1) tan" 1 (q£)] 
(13) 

where T is the Gamma function, r s is the spheres ra- 
dius, and £ gives the spatial extension of the aggregate. 
Incidentally, when D = 2, Eq. (JT3J yields S(q) = 
1 + 2(£/r s ) 2 (l + q 2 ^ 2 ) , which entails the sum of a 
Kerr contribution and a non-local one with Lorentzian 
lineshape. In the limit q£ << 1, Eq. (|13f) yields 



1 



D{D 



and the corresponding model for p reads as: 



1 - 



D(D 



7 ?7e r(£> + l)(^/r.) r 
P= 2kB~T 1 



(14) 



(15) 



From Eq. (|15fl it is readily seen that p spatially decays 
(when / = 0) as the modified Bcsscl function K$ of ar- 
gument £y/ D{D + 1)/6ct which depends on the fractal 
dimension. Recalling that a — r/wo and that the degree 
of nonlocality is the ratio between the spatial decay rate 
of the optically induced index perturbation and that of 
the self-trapped beams, our result implies that the de- 
gree of optical non-locality scales basically as the fractal 
dimension of the material. Noteworthy, Eqs. JSJ-^SJ) de- 
fine another well known model for non-local OSS, which 
applies in the case of nematic liquid crystals |2^| . 

In order to show that OSS exist also in the general 
case [Eq. <|13f) ]. with features directly linked to the fractal 
dimension D, we resort to numerical integration of Eqs. 
()7I8I) [or Eqs. (|10I8|) in the non-par axial regime] using 
finite difference discretization in a and Newton-Rapson 
iterations. To fix the ideas, we show results for the 
characteristic values of the following length scale ratios 
(between aggregate dimension £, particle radius r s , beam 
width w ): £/r s = 100 and K = £/w = 0.1. In Fig. □ 
(a) we show existence curves, i.e. the soliton normalized 
peak intensity against the normalized soliton width 
[max(u 2 ) vs. std of u 2 parametrized by /?], obtained for 
three values of D. As shown the features of OSS change 
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with fractal dimension D. This is even more clear from 
Fig-ffl(b), where we display the OSS normalized power 
(i.e. the norm Q = 2ir J u 2 rdr) and width as a function 
of the fractal dimension D (here we fix p = 2). In Fig. [21 
we show the effect of non-paraxiality. The features of 
OSS starts to exhibit significant deviations when the 
soliton width decreases significantly, and non-paraxial 
effects are no longer negligible. 

It is interesting to observe that for a fractal medium 
the equation for p in the configurational space includes 
fractional derivatives. Soft matter, not only provide a 
un-precedented framework to study nonlocality with a 
taylorable structure factor, but also opens the may to 
new mathematical models, as we will discuss in future 
publications. 

To prove that stable self-trapping can be achieved for 
input conditions that do not exactly match the OSS pro- 
file, we have run beam propagation simulations for the 
paraxial and non-paraxial models. In Fig. |21we show the 
spatial evolution of an input Gaussian TEMoo laser beam 
u(cr, £ = 0) = Acxp(— a 2 ), whose parameters do not ex- 
actly match the existence condition. The resulting longi- 
tudinal beam oscillations, which depend on beam power 
(for a fixed width) and degree of non-paraxiality, have 
been reported previously for other non-local solitons [sO] ■ 
They are connected with the fact that non-local solitons 
are absolutely stable, |2^| and can be related to excitation 
of "internal modes" of the soliton or, in the frame- 

work of the hyghly non-local approximation [i.e. when 
I(q) varies on a q-scale much broader than that of S(q)] 
to the existence of exact breathing solutions [13, HE E3 • 
Notably, in the latter regime, I(q) ~ 1(0) = P in Eq. 
0, with P the optical beam power. In this case, Eq. Q 
takes the form of the linear Schrodingcr equation for a 
quantum particle in a 2D potential well with shape dic- 
tated by particle autocorrelation g(r) and power P, and 
OSS reduce to the bound states, which can be found by 
standard techniques. Such results let us envisage a broad 
physical setting for the observation of deeply oscillating 
nonlocal solitons and show that their existence is not re- 
stricted to the paraxial regime, 

In summary we have shown that optical beams can 
be self-trapped in soft matter both in the paraxial and 
tightly focusing regimes, thus opening new perspectives 
for their applications in biological materials and as a 
mean to probe properties of condensed matter. For ex- 
ample, fixing the input beam waist and adjusting the in- 
coming power to find the appearance of a soliton allows 
one to directly measure the fractal dimension of the ag- 
gregates. Assessing the role of time dynamics and ther- 
mal contributions as well as the validity of the present 
approach in other condensed matter systems (e.g., su- 
percooled liquids, where optical trapping is unexplored 
to date) will be natural extensions of this work. 
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FIG. 1: (a)Normalized peak intensity max(u 2 ) Vs. the nor- 
malized beam width (i.e. the std of m 2 ), obtained from Eqs. 
flSt coupled to Eq. for various fractal dimensions D. 

(b) Normalized beam waist and power Vs. the fractal dimen- 
sion D for (3 = 2. 
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FIG. 2: As in figure soliton normalized peak intensity (a) 
and power (b) Vs. normalized waist for the paraxial (e = 0) 
and non paraxial (e = 0.1) model (D — 2.5, £ = 100r s , k = 
0.1). 
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